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ABSTRACT 

XRF 030723 is the first X-ray flash (XRF) to show in its optical light curve (LC) a bump 
that has been interpreted as the signature of a supernova (SN). After subtracting the afterglow 
component from the observed optical LC of the XRF counterpart, the properties of the putative 
SN are constrained by means of synthetic LCs of core-collapse SNe. For the redshift range z ~ 0.3 
- 1, all possible models require a rather small mass of synthesized 56 Ni, i.e., M( 56 Ni) ~ 0.01 - 
0.3 Mq. The models used to describe the energetic SNe Ic associated with gamma-ray bursts 
(SNe 1998bw and 2003dh) are too massive for the observed LC. If the relation between ejected 
56 Ni mass and total ejecta mass established from models of various Type Ic SNe also holds for 
the putative SN in XRF 030723, the ejecta mass is constrained to be ~ 1 - 3 Mq and the kinetic 
energy < 1 x 10 52 erg. This corresponds to a progenitor with 15M© < Mms ^5 25M©. The SN 
therefore appears to have properties intermediate between a normal SN Ic like SN 19941 and a 
more energetic object like SN 2002ap. 

Subject headings: gamma rays: bursts - supernovae: general - nucleosynthesis - X-rays: general 



I. INTRODUCTION 

X-ray flashes (XRFs) are intense transient 
bursts of X-rays, similar to gamma-ray bursts 
(GRBs) but characterized by a stronger X-ray 
than 7-ray fluence (Heise et al. 2001). It has been 
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speculated that XRFs are indeed linked to GRBs. 
Various scenarios have been proposed: XRFs may 
be GRBs viewed sufficiently off-axis (Yamazaki 
et al. 2003), or more massive explosions than those 
that make GRBs, leading to a smaller Lorenz fac- 
tor for the jet (Dermer, Chiang, & Mitman 2000). 
However, the real nature of the observed difference 
between GRBs and XRFs is not yet known. 

Both GRBs and XRFs are extragalactic sources. 
A number of GRBs have now been localized, and 
they are all cosmological. Interestingly, the three 
nearest long GRBs ever localized are all associ- 
ated with spectroscopically confirmed supernovae 
(SNe). Although the case of GRB 980425/SN 1998bw 
(Galama et al. 1998) was controversial, that of 
GRB 030329/SN 2003dh is established beyond 
doubt (e.g., Matheson et al. 2003. GRB 031203 at 
z = 0.1055 also appears to have hosted a SN pos- 
sibly similar to SN 1998bw (e.g., Thomsen et al. 
2004). 

All these SNe are of Type Ic (SNe Ic), and dis- 
play very high-velocity ejecta. SNe Ic are thought 



to originate from the collapse of the cores of mas- 
sive stars that have lost both their H and He en- 
velopes, exploding as C+O stars (e.g., Nomoto 
ct al. 1994). SNe 1998bw and 2003dh have been 
successfully modelled as highly energetic explo- 
sions (the spherical kinetic energy is ~ 30 — 50 
times that of a normal SN: Iwamoto et al. 1998; 
Mazzali et al. 2003), ejecting ~ 10 M of material 
and synthesizing ~ 0.5 M Q of 56 Ni, much more 
than in normal core-collapse SNe. The aspher- 
ical models require lower energy, but still 2—10 
times higher than a normal SN (Hofiich et al. 1999; 
Maeda et al. 2003). These SNe Ic and others, not 
known to be associated with GRBs but showing 
the spectroscopic signatures of a high explosion 
energy (> 10 52 erg), have been called hypernovae. 
Both the kinetic energy and the mass of 56 Ni of 
SNe Ic may positively correlate with progenitor 
mass (Nomoto et al. 2003). 

On the other hand, only a few XRFs have been 
accurately localized (to less than a few arcmin) 
thus far, and only one of these is known to be at 
low redshift: XRF 020903 at z = 0.251 (Soder- 
berg et al. 2003). In the case of XRF 030723, 
only upper and lower redshift limits could be de- 
termined, because the host galaxy could not be 
observed (Fynbo et al. 2004). The lower limit, 
z > 0.3, follows from the non-detection of the host, 
while the upper limit, z < 2.3, was derived from 
the lack of Lya absorption down to 4000 A. 

Fynbo et al. (2004) obtained optical photome- 
try and spectroscopy of XRF 030723. The i?-band 
light curve (LC) of the XRF counterpart showed a 
'bump', which may be the signature of a SN com- 
ponent. Other interpretations are possible (e.g., 
a two-jet model Huang et al. 2004), but Fynbo 
ct al. (2004) claimed that they could rule them 
out based on the SED evolution. Interpreting the 
bump as a SN sets further limits on the redshift 
(z ~ 0.3 - 1.0). Fynbo ct al. (2004) compared the 
bump with the LCs of different SNe Ic, concluding 
that the best match was given by the rest-frame 
B-band LC of SN 19941, at a redshift z = 0.6. 
SN 19941 was a normally-energetic SN Ic, with 
low-mass ejecta (Nomoto et al. 1994; Iwamoto 
ct al. 1994). However, the featureless spectrum of 
the bump obtained by Fynbo et al. (2004) makes 
them favour a broad-lined, highly energetic SN, 
like SN 1998bw or SN 2002ap. 



In this paper, we model the observed LC us- 
ing different SN Ic models, both normally- and 
hyper-energetic, and explore various redshifts to 
determine the possible range of the parameters of 
the SN associated with XRF 030723. 

2. Light Curve Models 

Assuming that the SN in XRF 030723 is of 
Type Ic, by analogy with the cases of GRB- 
associated SNe, we computed synthetic UVOIR 
bolometric LCs of exploded C+O star models. 
We used an LTE radiation hydrodynamical code 
and a gray 7-ray transfer code (Iwamoto et al. 
2000). The electron-scattering opacity was calcu- 
lated from the solution of the Saha equation, while 
the line opacity was fixed at 0.03 cm 2 g _1 , the 
value that was used to model SN 1998bw (Naka- 
mura et al. 2001). The total opacity is ~ 0.03 - 
0.1 cm 2 g _1 (compared to the range 0.05 - 0.15 
cm 2 g- 1 for SN 19941: Iwamoto et al. 1994). The 
accurate line opacity is unknown, if it is bigger, 
the timescale of LC peak becomes longer. 

We started from four different C+O star explo- 
sion models: C021, a normal SN Ic model devel- 
oped for SN 19941 (Iwamoto et al. 1994); CO100, 
an energetic model developed for SN 1997ef 
(Mazzali, Iwamoto, & Nomoto 2000); CO100/4, 
a scaled-down version of CO100 developed for 
SN 2002ap (Mazzali et al. 2002); and C0138, a 
model for the very energetic SN 1998bw (Naka- 
mura et al. 2001). The model parameters, i.e. 
ejected mass M e j, kinetic energy Ek, and progen- 
itor main-sequence mass Mms are summarized in 
Table 1. 

Given the uncertain redshift, we considered 
three values: z = 0.3, 0.6, and 0.8. The observed 
i?-band roughly corresponds to the rest-frame V- 
band if z = 0.3, the 5-band if z — 0.6, and the 
[/-band if z = 0.8. Because of time-dilation, the 
observed LC corresponds to narrower rest-frame 
LCs for increasing redshifts. On the other hand, 
in SNe Ic the LC peak is narrower for bluer bands. 
Coincidentally, these two effects roughly cancel 
out. 

We scaled the masses of the above C+O mod- 
els, conserving their explosion energies, and cal- 
culated synthetic LCs to find which scalings yield 
the best-fit models. Since the code computes bolo- 
metric LCs, the monochromatic LCs were esti- 
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mated using bolomctric corrections (BCs). As 
spectral and color information about the possible 
SN is very limited, we used the BCs of known SNe 
Ic as templates. These were calibrated in time 
with respect to the LC peak to match the appar- 
ently fast evolution of the SN. We used the BC 
template appropriate for the parameters of each 
model. For low-mass, low-energy models, we used 
the BCs of SN 19941 (assuming E(B - V) = 0.45; 
Richmond et al. 1996; B. Schmidt & R. Kirshncr 
1994, private communication). In other cases, the 
model parameters indicated energetic explosions, 
and then we adopted the BCs of a suitable hyper- 
nova and recalculated the model LC. In order to 
match the observed flux, we modified the mass of 
56 Ni synthesized for different redshifts and models. 
56 Ni powers the SN LC through the decay chain 

56 Ni ^ 56 Co ^ 56 Fe 

Model parameters affect the synthetic LC. The 
timescale of the LC peak, T pe ak, is determined by 
the combination of photon diffusion and ejecta ex- 
pansion and depends on Ek and M e j as T poa k « 
AK op t 1/2 M j 3/4 £ K ~ 1/4 (Arnctt 1982), where A 
represents the effect of the density structure and 
56 Ni distribution and K op t is the model-dependent 
average opacity. 

The parameters of our best-fit models are listed 
in Table 2, as are the template BCs used, and the 
computed synthetic LCs are compared to the ob- 
served LC in Figures 1, 2, and 3 for the three 
redshifts. The observed LC was derived from 
the Fynbo et al. (2004) R-band photometry by 
subtracting the Beuermann extrapolation of the 
power-law component of the afterglow. As ex- 
pected, for every value of i?K (i-e., for every model) 
there is a value of the ejected mass for which the 
synthetic LC reproduces the observations. Had 
we been able to model the spectra, we could have 
uniquely constrained the value of the kinetic en- 
ergy. However, the spectral information is limited. 
Fynbo et al. (2004) discussed that the SED of the 
bump is consistent with that of a SN similar to 
SN 1998bw, and that the only, low signal-to-noise 
spectrum is better consistent with a broad-lined 
SN like SN 1998bw or SN 2002ap than a narrow- 
lined one like SN 19941. We discuss these points 
later. 

Because the LC peak of the SN in XRF 030723 
is intrinsically narrow, the only models that give 
satisfactory results without major modifications 



to their mass values are models C021, which is 
scaled up somewhat, and CO100/4, which has 
to be scaled down (the top panel of Figure 4). 
All other models were originally developed for the 
high-mass HNe and must be scaled down in mass 
significantly to yield a narrow LC peak. 

The delayed LC-rise of the SN in XRF 030723, 
compared with the hyper-energetic SNe 1998bw 
and 2002ap (e.g., Galama et al. 1998; Yoshii et al. 
2003), cannot be reproduced by models with uni- 
form 56 Ni mixing. Instead, we had to assume that 
56 Ni is distributed only in the inner 10% of the 
ejected mass in our modelling. This increases the 
photon diffusion time, leading to a delayed onset 
of the LC. The total 56 Ni mass is ~ 0.012 - 0.015 
M Q if z = 0.3, - 0.07 - 0.12 M Q if z = 0.6, and 
~ 0.1 - 0.3 M Q if z = 0.8. These estimates would 
change if we could take into account the unknown 
if -corrections. For example, the if-correction of 
SNe la between the observed R- and rest-frame 
B-bands, Kbr, at z = 0.6 is ~ —0.6 (Kim, Goo- 
bar, & Perlmuttcr 1996). If a similar value was 
applied here (SNe Ic show some overall spectral 
similarity to SNe la), the 56 Ni mass for z = 0.6 
would decrease to ~ 0.04 M . In the case z — 0.3 
the relevant if -correction for SNe la, Kvr, varies 
between —0.2 and —0.5 near the peak, which sug- 
gests that the 56 Ni mass could decrease to < 0.01 
Mq. The if-correction for z = 0.8, Kjjr, is pos- 
sibly positive for SNe Ic because of the very small 
flux in the [/-band. This could make the 56 Ni mass 
larger than 0.25 M . 

3. Discussion 

All the models discussed above can reproduce 
the observed LC. Given the limited spectral in- 
formation, it is unfortunately not easy to select 
among them. However, based on their properties, 
we can at least attempt to narrow down the range 
of possibilities. 

In Figure 4, the basic parameters of the best-fit 
models are compared with the properties of four 
well-studied SNe Ic of various energies. The top 
panel is a plot of M j vs. Ek for the various mod- 
els and SNe. The bottom panel shows how M( 56 Ni) 
varies depending on the assumed redshift, all other 
model properties being the same. Redshift affects 
both the estimated SN luminosity and the BCs 
adopted in the LC modelling. Based on this plot, 
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we can select the models that are likely to be ap- 
plicable to the SN in XRF 030723, assuming that 
it is a SN Ic. The models for the observed Type 
Ic SNe are also plotted. 

First we discuss the low-redshift case (z = 0.3). 
In this case, the mass of 56 Ni is very small (< 
O.O15M ) for all models. This value is much 
smaller than even in a normal SNe Ic like SN 19941, 
suggesting that z — 0.3 is a significant underesti- 
mate of the real redshift if the observed bump is 
in fact a SN. However, if XRF 030723 really oc- 
curred at z = 0.3, then strong fall-back and/or a 
low-energy SN must be invoked to explain such a 
small 56 Ni mass. If fall-back is strong, much of the 
synthesized 56 Ni may be captured by the compact 
remnant, like in some Type II SNe which eject 
very little 56 Ni (e.g., SN 1997D, Turatto et al. 
1998; Zampieri et al. 2003). The other possibil- 
ity is that the explosion energy of the SN is too 
low to synthesize much 56 Ni. The Ek - M e j rela- 
tion of our best-fit models suggests that SNe with 
(M e j, Ek) as low as, e.g., (0.8, 0.5) would have 
the correct LC and may be expected to synthesize 
little 56 Ni as is required for z = 0.3. Such SNe 
must come from low-mass C+O stars, originating 
from main-sequence stars of < 15M Q . 

Next we consider the high-rcdshift case (z = 
0.8). This sets an upper limit for the ejected 56 Ni 
mass at M( 56 Ni) ~ 0.3 M Q . Combined with the 
tendency of actual SNe Ic to produce more 56 Ni for 
increasing M e j and Ek, this in turn sets an upper 
limit to M e j and Ek- The ejecta mass of the SN 
in XRF 030723 is probably less than ~ 6M© (i.e., 
Mms < 30 M Q ), and the kinetic energy less than 
~ 3 x 10 52 erg. These values are smaller than in the 
GRB-associated SNe. It is very unlikely that the 
SN can be as massive as SNe 1998bw (Iwamoto 
et al. 1998) or 2003dh (Mazzali et al. 2003) for 
two reasons: first, the kinetic energy would be ex- 
tremely high (~ 1 x 10 53 erg), and secondly, it 
would be strange for such a hyper-energetic ex- 
plosion only to synthesize ^0.2 M Q of 56 Ni. 

If the redshift z is <~ 0.6, the required 56 Ni 
mass is - 0.07 - 0.12M Q . Both the SN 19941- 
likc model *C021 (M oj - 1M Q , £ K ~lx 10 51 
erg, M( 56 Ni) ~ O.O7M ) and the SN 2002ap-like 
model *CO100/4 (M ej - 1.7M Q , £ K ~3x 10 51 
erg, M( 56 Ni) ~ 0.1M Q ) are consistent with the ob- 
served relations between M j, Ek, and M( 56 Ni). 
If the redshift is slightly higher (e.g., z ~ 0.7), a 



larger M( 56 Ni) is required, and a more massive 
and energetic model (M e j ~ 2 — 3M Q ; Ek ~ 4 - 
10 x 10 51 erg) is probably necessary. This makes 
the tentative SN more like SN 2002ap. Summa- 
rizing, our preferred estimates are for a redshift 
z ~ 0.5 — 0.7, consistent with the range favored by 
Fynbo et al. (2004), and SN properties between a 
normal SN and a low-energy, low-mass hypernova. 

These results are consistent with the colors of 
the putative SN on August 6 and 19, i.e., 14 
and 21 days after the XRF (Fynbo et al. 2004). 
We made a rough comparison between the XRF 
and the color evolution of other Type Ic SNe for 
the three redshifts, neglecting the unavailable K- 
corrections. The observed V — R ~ 0.1 — 0.25 and 

i — R ~ —0.4 0.2 on August 6 cannot be used 

to distinguish between the three redshifts and be- 
tween SNe 1998bw, 2002ap, and 19941. However, 
the U B non-detection on August 16, (B — R > 2 
and U - R > 2) seems to match SN 2002ap bet- 
ter because the U/UV flux is more depressed in 
SN 2002ap than in SNe 19941 or 1998bw. The 
detected K = 21.2 on Aug 14, however, suggests 
R — K ~ 3, which is too big to be consistent with 
the photometry of SNe 1998bw, 2002ap, and 19941 
(but see Fynbo et al. 2004). Finally, the observed 
i — R ~ —1.1 on August 16 favors z ~ 0.6 and 
a comparison with SNe 19941 and 2002ap. Sup- 
port for an SN 2002ap-like object also comes from 
the apparently featureless spectrum of the bump 
(Fynbo et al. 2004). SN 2002ap has a broad-lined 
spectrum, much more similar to SN 1998bw than 
to SN 19941. 

This is the first case of an XRF probably as- 
sociated with a SN (Watson et al. 2004). Inter- 
estingly, the properties of the SN appear to lie 
at the low end of the distribution of HN prop- 
erties, or perhaps even to be similar to those of 
normal SNe Ic (Figure 4). This may be the con- 
sequence of an unfavorable orientation which pre- 
vented us from seeing the GRB and the most en- 
ergetic part of the ejecta, but it also may indi- 
cate a real difference between XRFs and GRBs 
on the one hand and the accompanying SNe on 
the other. If the progenitor of the SN/XRF is 
really intermediate between those of SN 2002ap 
and 19941, i.e., 15M Q < M MS < 25M Q , a black 
hole may not be formed (e.g., if Mms ^ 20M Q 
as in Fryer & Kalogera 2001), and the central en- 
gine may be a neutron star like in the "magnetar" 
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model (Nakamura et al. 1998), instead of the black 
hole of the "collapsar" model for long GRBs (Mac- 
Fadyen & Woosley 1999). Accurate spectral infor- 
mation, which would have allowed us to constrain 
the properties of the SN much more tightly, is un- 
fortunately not available. There is however one 
case where a normal, or possibly SN 2002ap-likc 
SN Ic was claimed to be associated with a GRB 
(GRB 021211/SN 20021t, Delia Valle et al. 2003). 
Clearly, more data on the SN - GRB/XRF con- 
nection are necessary before we can understand 
the full extent of the relation between these phe- 
nomena. 
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Table 1 
Original SN/HN Models 



Model 


Supernova 


M MS /M 


M oj /M 


E K 








(10 51 ergs) 


C021 


SN 19941 


~ 15 


0.88 


1 


CO100/4 


SN 2002ap 


20 - 25 


2.4 


4.8 


CO100 


SN 1997ef 


35 - 40 


9.6 


20 


C0138 


SN 1998bw 


~ 40 


11 


50 



Table 2 
Best-fit Models 



Model a 


M ej /M Q 


E K 




M( 56 Ni)/M 




Bolomctric Correction 






(10 51 ergs) 


z = 0.3 


z = 0.6 


z = 0.8 


Template b 


*C021 


1 


1 


0.012 


0.07 


0.10 


SN 19941 


*CO100/4 


1.7 


4.8 


0.014 


0.11 


0.27 


SN 2002ap 


*CO100 


3.2 


20 


0.015 


0.12 


0.30 


SN 2002ap 


*CO!38 


6.4 


50 


0.014 


0.09 


0.12 


SN 1998bw 



a The subscript * is used to discriminate each modified model from the original model. 

b The BCs for SNe 1998bw and 2002ap arc scaled up in time by a factor of 1.5 and 1.2, respectively, with 
respect to the peak epoch. 



6 




days since XRF (rest frame) 



Fig. 1. — Comparison between the best-fit model LCs (solid line: *C021; long-dashed line: *CO100/4; 
short-dashed line: *CO100; dotted line: *C0138) and the observed LC for z=0.3 (filled squares; Fynbo et al. 
2004). 
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Fig. 2. — Same as Fig. 1, but for z—0.6. 
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Fig. 3. — Same as Fig. 1, but for z=0.8. 
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Fig. 4. — Top panel: M e j vs. Ek for the best-fit models (squares) and for four well-studied SNe Ic (filled 
circles). Bottom panel: M e j vs. M( 56 Ni) for the best-fit models (squares) at z = 0.8 (top), 0.6 (middle), and 
0.3 (bottom), respectively, and for the SNe Ic (filled circles). 



10 



